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In the mirror world hypothesis the mirror baryonic component emerges as a possible dark 
matter candidate. An immediate question arises: how the mirror baryons behave and what 
are the differences from the more familiar dark matter candidates as e.g. cold dark matter? 
In this paper we answer quantitatively to this question. First we discuss the dependence of 
the relevant scales for the structure formation (Jeans and Silk scales) on the two macroscopic 
l/^ ' parameters necessary to define the model: the temperature of the mirror plasma (limited by 

O^l , the Big Bang Nucleosynthesis) and the amount of mirror baryonic matter. Then we perform 

a complete quantitative calculation of the implications of mirror dark matter on the cosmic 
microwave background and large scale structure power spectrum. Finally, confronting with 
Xjr^ ' the present observational data, we obtain some bounds on the mirror parameter space. 
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cn ! 1 Introduction 
O ^ 

r~| ', The idea that there may exist a hidden mirror sector of particles and interactions with exactly the 
Oih| same properties as that of our visible world was suggested long time ago . The basic concept 
Q ■ is to have a theory given by the product G x G' of two identical gauge factors with the identical 
particle contents, which could naturally emerge e.g. in the context of Eg x Eg superstring. (From 
now on the "primed" (') fields and quantities stand for the mirror (M) sector to distinguish from 
the ones belonging to the observable or ordinary (O) world.) In particular, one can consider a 
- minimal symmetry GsM x G'o^m where Gsm = SU{3) x SU{2) x U{1) stands for the standard 
r> I model of observable particles: three families of quarks and leptons qi, u^, d?; k, ef (i = 1,2,3) 
' and the Higgs doublet (j), while Ggp^^ = [SU{3) x SU{2) x U{1)]' is its mirror gauge counterpart 
with analogous particle content: fermions q'^, u'f, d'f; l[, ef and the Higgs </>'. The M-particles 
are singlets of Gsm and vice versa, the 0-particles are singlets of GgM- More generally, one can 
have in mind the grand unified extensions like SU{5) x SU{5y, 50(10) x S'O(IO)', etc. Besides 
the gravity, two sectors could communicate by other means In particular, ordinary and 
mirror neutral particles could mix: e.g. photons |H] or neutrinos [1], two sectors could have a 
common gauge group of flavour jS] or common Peccei-Quinn symmetry [Hj. 

A discrete symmetry G <^ G' interchanging corresponding fields of G and G' , so called 
mirror parity, guarantees that two particle sectors are described by identical Lagrangians, with 
all coupling constants (gauge, Yukawa, Higgs) having the same pattern. As a consequence the 
two sectors should have the same microphysics. In the particular case in which G sector is 
left-handed and G' sector is right-handed, this discrete symmetry can be interpreted as the true 
parity 
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If the mirror sector exists, then the Universe along with the ordinary photons, neutrinos, 
baryons, etc. should contain their mirror partners. One would naively think that the O- and M- 
sectors, having the same microphysics, should have also the same cosmology.^ Then one would 
expect that M-particles are present in the same amount in the Universe as 0-ones, which would 
be in immediate conflict with the Big Bang Nucleosynthesis (BBN) bounds on the effective 
number of extra light neutrinos AN^: the mirror photons, electrons and neutrinos would give a 
contribution equivalent to AA^j^ ~ 6.14. However two sectors may have different initial condi- 
tions. In particular, after inflation, the two sectors could be reheated at different temperatures, 

7^ Tr, which can be achieved in certain inflationary models P^-jllj. If the O- and M-particles 
interact weakly enough (which condition is automatically fulfilled if the two worlds communi- 
cate only via the gravity), the two systems do not come into thermal equilibrium at later epoch 
and they evolve independently, maintaining approximately constant the ratio x = T' /T among 
their temperatures. The BBN constraints are satisfied if x is sufficiently low. Namely, the most 
conservative bound ANy < 1 [12 implies that x < 0.64 jllj . 

The difference of the temperatures T' < T breaks the symmetry between the cosmological 
properties of O- and M-sectors. Namely, the number density of M-photons is much smaller than 
that of 0-ones, n'^/n^ = x^ <C 1. It is important to understand, whether the mirror baryon 
density could be comparable or even larger than the ordinary one Ub, in which case they 
could constitute dark matter of the Universe, or at least its significant fraction. 

The cosmological evolution of the Mirror Universe was studied in ref. 11 and all the key 
epochs (baryogenesis, nucleosynthesis, recombination etc.) were analyzed in details. It was 
shown, in particular, that in the context of the GUT or electroweak baryogenesis scenarios the 
condition T' < T yields that ry^ > ry^,, where % = Uh/n^ and r/^ = 11^11'^ respectively are the 
baryon asymmetries in O- and M-sectors |llj . However, t/j,///;, > 1 is not yet enough for having 
f3 = nyrih > 1. Since P = x^^rj^-qh), the latter requires much stronger condition r/^/r/b > x~^. As 
it was demonstrated in ref. this condition can be satisfied for a certain range of parameters, 
with the values of x which are not too low, x > 0.01 or so. However, more appealing situation 
emerges in a leptogenesis scenario due to particle exchange between the ordinary and mirror 
sectors suggested in ref. which predicts /? > 1 but also implies that /3 < 10 or so, and thus 
can explain naturally the near coincidence between the visible (0-baryon) matter density 0,^ 
and the dark (M-baryon) matter density fij^ in a rather natural way [2]. 

If > niirror baryons emerge as a possible dark matter candidate; they can contribute 
the dark matter of the Universe along with the cold dark matter or even constitute a dominant 
dark matter component. An immediate question arises: how the mirror baryon dark matter 
(MBDM) behaves and what are the differences from the more familiar dark matter candidate 
as the cold dark matter (CDM)? 

The peculiar properties of mirror dark matter were discussed qualitatively in and this 
analysis was confirmed and extended in refs. |14[ I15j. In this paper we complete this program 
giving a complete quantitative presentation of the implications of the MBDM for the cosmolog- 
ical large scale structure (LSS) and the cosmic microwave background (CMB) anisotropies. 

The plan of the paper is as follows. In the next section we analyze the relevant length scales 
for structure formation in the mirror photon-baryonic sector. In section 3 we describe the effect 
of the evolution of perturbations in linear regime and compute the power spectra for the LSS 
and CMB for various values of the two parameters x and /?. The main differences with respect 
to a standard CDM scenario are discussed. Finally, our main conclusions are summarized in 
section 4. 

^More generically, M-parity could be spontaneously broken and the weak scales {(f>) and {((>') could be different, 
which would lead to somewhat different particle physics in the mirror sector |S]. 
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2 Relevant length scales 



Mirror matter may seem a tremendously complicated dark matter candidate. However, from the 
point of view of structure formation, it can be described relatively simply. The microphysics of 
the mirror sector is in fact well defined, being identical to that of our sector. All the differences 
with respect to the ordinary world can be described in terms of two macroscopic parameters 
which are the only free parameters in the model: 

T' Vt' 

^^j^ (1) 

where T (T') is the ordinary (mirror) photon temperature in the present Universe,^ and Qb (fi^,) 
is the ordinary (mirror) baryon density fraction. In this section, we discuss the dependence of 
the length scales relevant for structure formation from these parameters. 

In the most general context, the present energy density contains relativistic (radiation) 
component Qr, non-relativistic (matter) component and the vacuum energy density Q\ 
(cosmological term). The present observational data indicate that the Universe is almost flat, 
i.e. = + + ~ 1 in a perfect accordance with the inflationary paradigm, with 
= 0.2 — 0.3 and the rest of the energy density is due to the cosmological term. In the 
context of our model, the relativistic fraction is represented by the ordinary and mirror photons 
and neutrinos, Q^h'^ = 4.2 x 10~^(1 + x^), where contribution of the M-species is negligible due 
the BBN constraint <^ 1. As for the non-relativistic component, it contains the 0-baryon 
fraction Qf, and the M-baryon fraction 0,'^ = Pil-b, while the other types of dark matter, e.g. the 
CDM, could also present and so in general = fif, + + ficDM-^ In the following we use 
the central values by WMAP LUb = Vlbh? ~ 0.023 and oJm = ^mh'^ ~ 0.135 |19j . and consider 
scenarios with /3 = 1 5 where the limiting case /? ~ 5 corresponds to the case when the dark 
matter is entirely due to MBDM (Ocdm = 0). 

The important moments for the structure formation are related to the matter-radiation equal- 
ity (MRE) and to the plasma recombination and matter-radiation decoupling (MRD) epochs. 
The MRE occurs at the redshift: 

3240 f LOm \ 

' + 17, ^ 1 + [0.135) ^ > 

where we denote = ^mh'^- Therefore, for x ^ 1 it is not altered by the additional relativistic 
component of the M-sector. 

The matter radiation decoupling takes place only after the most of electrons and protons 
recombine into neutral hydrogen and the free electron number density strongly diminishes, so 
that the photon scattering rate drops below the Hubble expansion rate. In the ordinary Universe 
the MRD takes place in the matter domination period, at the temperature Tjec — 0.26 eV which 
corresponds to redshift 1 -|- ^^ec = Tdec/Ttoday ^ HOO. 

The MRD temperature in the M-sector T^^^ can be calculated following the same lines as in 
the ordinary one Due to the fact that in either case the photon decoupling occurs when the 
exponential factor in the Saha equation becomes very small, we have T^^^ ~ Tdec, up to small 
logarithmic corrections related to r]' , different from rj. Hence 

/ -1. X 1100 

1 + ^dec -X (1 + Zdec) ^ (3) 

X 



^The ratio of the temperatures between two sectors is nearly constant during the evolution of the Universe. 
In general, one has T'/T = x[gs(T) / g's{T')]^^^ , where the factors gs and g'^ accounting for the degrees of freedom 
of the two sectors can be different from each other, see ref. for details. 

^In the context of supersymmetry, the CDM component could exist in the form of the lightest supersymmetric 
particle (LSP). It is interesting to remark that the mass fractions of the ordinary and mirror LSP are related 
as il'i^sp — xHlsp- In addition, an HDM component fli, could be due to neutrinos with mass < 0.3 eV. The 
contribution of the mirror neutrinos scales as Q'^ = x'^Q,„ and thus it is irrelevant. 
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so that the MRD in the M-sector occurs earher than in the ordinary one. Moreover, for a value 
X = Xeq, where: 

X = 34 f (4) 

=4 1 + 

the mirror photon decouphng epoch coincides with the MRE epoch. This critical value plays 
an important role in our further considerations. Namely, for x < Xeq the mirror photons would 
decouple yet during the radiation dominated period. 

Let us discuss now the relevant scales for evolution of perturbations in the MBDM. The 
relevant scale for gravitational instabilities is Jeans length, defined as the minimum scale at 
which, in the matter dominated epoch, sub-horizon sized perturbations start to grow. The 
mirror Jeans scale is given by: 

X'jiz)^v',{z)i^/Gp{z)y/'{l + z) (5) 

where p{z) is the matter density at a given redshift z, v'g{z) is the sound speed in the M-plasma 
and the (1 + z) factor translates the physical scale at the time of redshift (1 + z) to the present 
scale. We remark that the M-plasma contains more baryons and less photons than the ordinary 
one, pj, = Ppb and p'^ = x'^p^, and thus the sound speed can have a quite different behaviour. 
We have: 




-1/2 



1 

73 



. , 3/3u;6 _4 fl + Zeq 

1 -I- (1 -I- X ^ ' 



4 \ 1 + Z 



-1/2 



(6) 



where ujf, = O^/i^, quite in contrast with the ordinary world, where Vs ~ c/\/3 practically till the 
photon decoupling. The M-baryon Jeans length reaches the maximal value at z = z'^^^, where 
it is given by ^: 

100 x^/^ 

After decoupling, eq. © does not hold anymore and the Jeans scale decreases to very low 
values, due to the fact that the pressure supplied by the relativistic component of the mirror 
plasma disappears. 

Density perturbations in MBDM on scales A > Aj which enter the horizon at z ~ z^q 
undergo uninterrupted linear growth immediately after z^q. Perturbations on scales A < Aj j^^, 
start instead to oscillate immediately after they enter the horizon, thus delaying their growth 
till the epoch of M-photon decoupling. 

Finally, we turn our attention to dissipative processes which can modify the purely gravita- 
tional evolution of perturbations. As occurs for perturbations in the 0-baryonic sector, also the 
M-baryon density fluctuations should undergo the strong collisional damping around the time 
of M-recombination. The photon diffusion from the overdense to underdense regions induces a 
dragging of charged particles and washes out the perturbations at scales smaller than the mirror 
Silk scale, A^. The behavior of A^. as a function of the parameter x and /3 is given by 

where f{x) = X^l"^ for X > Xeq, and /(x) = (x/Xeq)^/^Xeq'^ for X < Xeq. 

The impact of such scales on the evolution of density perturbations will be discussed in the 
next section. 



*We cissumed fix * S> 1 so that the constant term in eq. ® can be neglected. 
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3 Evolution of perturbations 



We clearly understand from the previous discussion that MBDM has peculiar features which 
can leave a characteristic imprint in the large scale structure of the Universe. 

First, perturbations in MBDM on scales A < Aj experience an oscillatory regime. The 
MBDM oscillations transmitted via gravity to the ordinary baryons, could cause observable 
anomalies in LSS power spectrum and in the CMB angular power spectrum. 

Second, for x > Xgq, the growth of perturbations on scales A < Aj ^g^, does not start at z^q 
but is delayed till M-photon decoupling. If MBDM is the dominant dark matter component, 
one expect to observe less structures on these scales than in standard CDM scenario. 

Finally, the density perturbation scales which can run the linear growth after the MRE 
epoch are limited by the length X'g. To some extent, the cutoff effect is analogous to the free 
streaming damping in the case of warm dark matter (WDM) . 

In order to make quantitative predictions we computed numerically the evolution of adiabatic 
perturbations in a Universe in which is present a significant fraction of mirror dark matter. More 
precisely, following the approach described in |17j . we solved numerically in a synchronous gauge 
the linear evolution equations for perturbations of all matter components: ordinary baryons, 
photons, neutrinos, their mirror analogues, and cold dark matter. In fact, with respect to the 
standard case, the full set of equations was doubled in order to properly take into account the 
evolution of the mirror photon-baryon system. The decoupling in ordinary and mirror plasma 
was followed numerically as prescribed in ^Ji- All computations were made assuming a flat 
space-time geometry (JIq = 1; i-e. ^^a = 1 — ^m)- In order to compare our predictions with the 
"standard" CDM results, we have chosen a "reference cosmological model" with scalar adiabatic 
perturbations and no massive neutrinos with the following set of parameters |16| : 



a;fe = 0.023, = 0.133 (17™ = 0.25), J7a = 0.75, 



0.97, h = 0.73. 



(9) 



We have included in this model the mirror sector and studied the effects of MBDM as a function 
of the parameters x and /3. For the sake of comparison, in all calculations the total amount of 
matter = ^^CDM+^^b+^b maintained constant, 0,.^ = 0.25. Mirror baryons contribution 
is thus always increased at the expenses of diminishing the CDM contribution. Hence, the case 
when dark matter is entirely represented by the MBDM, with no CDM component, corresponds 

to P = UJrn/^^b - 5. 
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Figure 1: Evolution of perturbations at the scale fc = 0.1 Mpc^^ in the case when dark matter is 
entirely due to mirror baryons (ficDM = 0). Dot-dashed and dotted lines correspond to ordinary baryons 
and photons, while long dashed and dashed lines are for mirror baryons and photons. All cosmological 
parameters are taken as in Left panel (a) corresponds to the case x = 0.6 and the right panel (b) to 
the case x — 0.3. 
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Figure 2: The same as in Fig. ^ but in the case when the MBDM is a sub-dominant dark matter 
component with (3 — 2, i.e. — 0.09 and ficDM = 0.12. 

To understand the impact of mirror matter on structure formation it is useful to look at Fig. 
1 and Fig. 2 where, for a selected wavenumber k = 27r/A and for selected values of x and j3, we 
show the evolution of the density contrast in the various components, 6i = 5pi/pi, as a function 
of the scale factor a. In Fig. 1 we show the situation when dark matter is entirely due to mirror 
baryons, while in Fig. 2 mirror baryons represent a subdominant component. The panels (a) 
of the two figures are obtained by assuming large values of the mirror to ordinary temperature 
ratio (i.e. x = 0.6), while the panels (b) correspond to the value x = 0.3 which, for the chosen 
cosmological parameters, implies that M-photon decoupling approximately coincides with MRE 
epoch (c.f. eq. Q)). 

As long as the perturbation scale is larger than the horizon, the various components grow 
at the same rate {6i oc a^). The situation drastically changes when the perturbations enter the 
horizon (around a/ao ~ 10^^). At this point, baryons and photons, for each sector separately, 
become causally connected and behave as a single fluid. In other words, we have two fluids, 
the ordinary baryon-photon and the mirror baryon-photon ones. The sub-horizon evolution of 
these fluids depends on the value of the Jeans lengths. If the Jeans length is larger than the 
perturbation scale then the photo-baryon fluid starts to oscillate. This is always the case (before 
the decoupling) for ordinary baryons and photons. This is not always true in the mirror sector. 
By comparing panels (a) and (b) of these two figures we see in fact that for large values of x the 
M-photons and baryons oscillate, while for smaller x values perturbations undergo uninterrupted 
growth even before M-photon decoupling (which occurs around a/ao ~ lO^^x^^). 

Oscillations in the mirror sector (when present), transmitted via gravity to the ordinary 
baryons, would produce observable anomalies in LSS power spectrum and in the CMB anisotropy 
spectrum. By comparing the late time perturbation evolution in Fig. la and Fig. 2a, we can 
understand that the efficiency of this process depends on the amount of mirror matter present 
in the Universe. After decoupling, in fact, baryons, which are no longer supported by photon 
pressure, rapidly fall in the potential wells created by the dominant dark matter component. 
If M-baryons dominate the dark matter budget, this leads asymptotically to 5b = (5^ (see late 
time evolution in Fig. la). This means that the baryonic structure power spectrum re- write 
the M-baryonic power spectrum, which is suppressed at small wavelength due to Silk damping 
and is eventually modulated (if x is not too small) as a results of acoustic oscillation. If instead 
CDM is the dominant dark matter component, we have asymptotically Sb = = Sqdm (see 
late time evolution in fig. 2a), which means that both M-baryonic and 0-baryonic structures will 
follow the "standard" CDM power spectrum. 

The dependence of the LSS power spectrum on the parameters x and (3 is shown explicitly 
in Fig. 3. In the upper panel we assume that the dark matter is entirely due to mirror baryons 
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and we consider variations of the x parameter. For large x values, as a result of the oscillations 
in MBDM perturbation evolution, one observes oscillations in the LSS power spectrum. The 
position of these oscillations depends on x, as can be easily understood. The smaller is x the 
smaller is the mirror Jeans scale at decoupling and thus the smaller are the perturbations scales 
which undergo acoustic oscillations. Superimposed to oscillations one clearly see the cut-off in 
the power spectrum due to the Silk damping. We remark that the Silk scale also depends on 
X (and f3), as it is described by eq. (jHJ. As a consequence, the cut off in the power spectrum 
moves to smaller wavelength when we decrease the x parameter. 

In the lower panel of Fig. 3 we can appreciate the role of the parameter (3. One clearly sees 
that, as expected, the smallest is the amount of mirror baryons the less evident are the features 
in the LSS power spectrum. Interestingly, for large x values, one observes a relevant effect even 
for relatively small amount of mirror baryons; even for /? = 1 the effects are quite noticeable. 
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Figure 3: LSS power spectrum in the linear regime for different values of x and uj[^ = ft[,h^, as compared 
with a standard model prediction (solid line) . In order to remove the dependences of units on the Hubble 
constant, we plot on the x-axis the wave number in units of h and on the y-axis the power spectrum in 
units of h~^. All parameters are taken as in Q. We also show the binned data of 2dF observations 1.18.1. 
Top panel. Models where dark matter is entirely due to MBDM (no CDM, i.e. /3 = 5) for different values 
X = 0.3, 0.5, 0.7. Bottom panel. Models with mixed CDM+MBDM components, /? = 1, 2, 3, 4 for a value 
of X = 0.7. 



Finally, in Fig. 4 we show the effects of MBDM on the angular spectrum of the CMB 
anisotropy. The predicted spectrum is quite strongly dependent on the value of x (see upper 
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panel), and it becomes practically indistinguishable from the CDM case for x < Xcq — 0.3, in 
this case the MBDM behaves just as the CDM at the scales relevant for the CMB oscillations. 
However, the effects on the CMB spectrum rather weakly depend on the fraction of mirror 
baryons (see lower panel for different values of In other words, the CMB anisotropy spectrum 
is mainly sensitive to amount of extra-radiation in the Universe due to the the mirror sector 
(which is fixed by x, being $7', oc x^) but it can hardly distinguish between MBDM and CDM. 
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Figure 4: The CMB angular power spectrum for different values of x and oj^, compared with a standard 
model (solid line). We also show the WMAP ^H] and ACBAR 20;, data. The choices of the parameters 
for both top panel and bottom panel exactly correspond to those of Fig. 3. 

Our predictions can be compared with the observational data in order to obtain bounds on 
the possible existence of the mirror sector. To give a visual impression of the present situation 
we show in the Fig. 3 the 2dF binned data and in Fig. 4 the WMAP ^ and ACBAR ^ 
data. 

To extract information from the experimental data one clearly needs a detailed statistical 
analysis. However, some general conclusions can be obtained in a rather straightforward way: 

(i) The assumption that DM is entirely due to mirror baryons is evidently not compatible 
with present LSS data unless the value of x is enough small: x < Xgq ~ 0.35. 

(ii) Very high values of x, x > 0.5, can be excluded even for a relatively small amount of 
mirror baryons. E.g. for x = 0.7, one has relevant effects on LSS and CMB power spectrum 
down to values of M-baryon density of the order $7^ ~ On the other hand, intermediate 



8 



values X = 0.3 — 0.5 can be still allowed if the MBDM is a subdominant component of dark 
matter. 

(iii) For small values of x, say x < 0.3, neither the linear LSS power spectrum nor the CMB 
angular power spectrum can distinguish the MBDM from the CDM. In this case, in fact, the 
Jeans length Aj and the Silk length A^, which mark region of the spectrum below which one 
sees the effects of mirror baryons, decrease to very low values, which undergo non linear growth 
from relatively large red-shift. 

4 Conclusions 

The idea of a mirror sector of particles and interactions has attracted a significant interest over 
last years. The concept of mirror world could have interesting implications for the following 
problems: detection of Machos via gravitational microlensing |H1 |^, search ordinary-mirror 
star binaries |221 and mirror planetary objects [23]) implications of mirror matter for galaxy 
halos 121] , the role of the sterile neutrinos in neutrino physics jlj |23] and their implications for 
ultra- high energy cosmic rays IQ. and gamma ray bursts |2S]) implications for the flavor and 
CP violation and axion physics 0101) etc. In particular, the effects based on the kinetic mixing 
among the ordinary and mirror photons [HI EJj can provide interesting possibilities for revealing 
the mirror matter in positronium decays [2H1 and in dark matter detectors .29_. 

In this paper, we have studied the cosmological implications of the mirror matter as dark 
matter. More precisely, we have quantitatively discussed the effect of the MBDM on the evolu- 
tion of density perturbations in the linear regime as function of the two free parameters in the 
model: the temperature of the mirror plasma (limited by BBN to a; = T'/T < 0.64) and the 
amount of mirror baryonic matter /? = QyQiy. We summarize here the main conclusions. 

The concept of mirror baryons, as possible candidate for dark matter, introduces two new 
scales in the structure formation scenario: the Jeans scale of the mirror photo-baryon fluid, Aj, 
and the Silk damping scale of mirror baryons, Ag. Due to the pressure support of mirror photons, 
perturbations in the mirror fluid on scales smaller than Aj cannot grow before mirror photon 
decoupling. If decoupling occurs after the MRE epoch (i.e. if x > x^q — 0.35), and if mirror 
baryons are a relevant dark matter component, one expect then to see less structures on these 
scales with respect to the standard CDM scenario. In addition, mirror baryon perturbations on 
scales A < Aj go through an oscillatory regime. This could produce, via gravity, observable 
distortions in the CMB and LSS power spectrum. Finally, the mirror Silk scale Ag introduces 
a cut-off for the perturbation scales which can run the linear growth after matter-radiation 
equality. 

All these effects were taken into account numerically in the Fortran code which we used to 
follow the perturbation evolution in presence of the mirror sector and to compute the LSS and 
CMB power spectra. The results are shown in Figs. 3 and 4 for various choices of parameters 
X and p. From a comparison with present observational data, one is able to conclude that the 
existence of a mirror sector with a relatively high temperature, T' > 0.4 T and a high mirror 
baryon density Jl^ > few can be already excluded. This confirms previous estimates of refs. 
im I14j and improves the BBN bound on the mirror sector parameter space. More details can 
be found in ref. jJH]. 

In order to further reduce the allowed parameter space, it would be clearly important to 
extend the analysis to the non-linear regime. Many relevant questions, related to the dynamics 
in this regime, can be formulated: what implications could have the mirror Silk cutoff for scales 
which already went non-linear, e.g. for scanning the early Universe via studying reionization, 
quasars, Lyman-a forest, etc.? Can mirror baryons, being a dissipative dark matter component, 
provide extended triaxial halos instead of being clumped into the galaxy as usual baryons? 
How the star formation mechanism proceeds in the M-sector where the temperature/density 
conditions and chemical contents are much different from the ordinary ones? Could the mirror 
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protogalaxy at a certain moment before disk formation become a collisionless system of the 
mirror stars and thus maintain a typical elhptical structure? What is the initial mass function 
of mirror stars, and how many and how heavy mirror stars do we expect as Machos in the 
galactic halo? 

Many other questions can be formulated and various new data are needed to discriminate 
different cosmological settings. Our numerical analysis, which describes the mirror world in the 
linear regime, sets the starting point for answering these questions. 
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